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CL Abstract. We introduce a new method for estimating the growth 

of various quantities arising in dynamical systems. We apply our 
method to polygonal billiards on surfaces of constant curvature. 
For instance, we obtain power bounds of degree two plus epsilon 
for billiard orbits between almost all pairs of points in a planar 
polygon. 
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Introduction and overview 



Complexity of a dynamical system is measured with respect to a 
coding of its orbits. The coding, in turn, is determined by partitioning 
the phase space of the system into elementary pieces. For dynami- 
i/-) ' cal systems with singularities, such as polygonal billiards, connected 

components in the complement to the singular set yield a natural par- 
tition. Convexity of its atoms with respect to the geodesic structure in 
the phase space imposed by geometric optics, is crucial in the study of 
O \ billiard complexity [7J. 

In the present study, as well as in [TJ, P is a geodesic polygon in 
a surface of constant curvature. Let, for concreteness, P be a planar 
polygon. We denote by fp(n) the number of words of length n gener- 
^ ' ated by coding billiard orbits by visited domains of regularity. When 

^ ■ P is simply connected, this coincides with the coding by sides in P. It 

is known that fp{n) is subexponential in n [31 E], and for general P 
no better bound is known. If P is a rational polygon (i. e., its angles 
are commensurable with 7r [3]), fp(n) = 0(n 3 ) (H [7J. The current 
conjecture is that for any planar polygon fp(n) = 0(n d ) [5]. 

In order to advance the understanding of billiard complexity, we 
introduce the notion of partial complexities. Let \1/ be the phase space, 
and let V be the defining partition. Iterating the dynamics we obtain 
an increasing tower V n of partitions; the full complexity is f(n) = 
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\V n \. If R C $, let V n (R) be the induced tower of its partitions. The 
partial complexity based on R is fain) = \V n (R)\. Particular partial 
complexities have been studied earlier. For instance, in [8] we obtained 
polynomial bounds on direction complexity, which is one of the partial 
complexities investigated here. 

In this work we introduce a new general approach to estimating par- 
tial complexities. The setting is as follows. There is a family of subsets 
Rq foliating the phase space. Let fe(n) be the partial complexity with 
base Rq. Let ge(n) be the counting function for singular billiard or- 
bits starting from Rq. Under appropriate assumptions, fQ(n) and gQ(n) 
have the same growth, as n — > oo. See section HJ 

Let 9 G 6, the parameter space. Suppose that we bound the average 
counting function G(n) = J e gg(n). Tchebysheff inequality and the 
zero-one law yield bounds for individual ge(n) valid for almost all 9 G 0. 
See section |2j Combined with preceding remarks, these yield estimates 
on partial complexities for almost all values of the parameter. 

This is the general scheme for our approach to partial complexities. 
This work implements this scheme for polygonal billiards. We will now 
describe the contents of the paper in more detail. 

In section[T]we investigate counting functions and their averages. We 
establish the relevant framework in sufficient generality, with the view 
towards a broad range of geometric-dynamic applications. The main 
results are Propositions [T] and [2] respectively. These yield geometric 
formulas for averages of counting functions which are valid under mild 
assumptions of transversality type. 

Section [2] is analytic, and also quite general. The setting is as follows. 
There is a family of positive functions, ge(p), of positive argument 
{p G N and p G M+ in the discrete and continuous cases respectively), 
depending on parameter 9 G O. Set G{p) = f e ge(p)d9. From upper 
bounds on G(p) we derive estimates on individual ge(p); they are valid 
for almost all 9. Precise formulations depend on the details of the 
situation. See Propositions [3] and HI 

Section [3] sets the stage for applications to billiard dynamics. Our 
billiard table is a geodesic polygon, P, in a simply connected surface 
of constant curvature. There are two versions of billiard dynamics: 
the billiard flow and the billiard map. In our discussion of partial 
complexities, it is convenient to treat them separately. Accordingly, 
section [3] consists of several subsections; each subsection deals with a 
particular partial complexity for a particular kind of billiard dynamics. 

We use two geometric parameters for partial complexities: the direc- 
tion and the position. The direction complexity tells us how the set of 
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phase points starting in the same direction splits after bouncing off of 
the sides of P. The direction complexity is defined for planar polygons. 
The position complexity tells us about the splitting of beams of billiard 
orbits emanating from a point of P. It is defined in all cases. 

In each of the subsections of section [3] we define a counting function 
and check the assumptions of section [IJ then we evaluate the integral 
over the parameter space, i. e., we compute the average counting func- 
tions. It turns out that they have geometric meanings. Here is a sample 
of results from sectional Let Gp{l) be the average position counting 
function for the billiard flow in a geodesic polygon P. For planar poly- 
gons we have Gp(l) = Cq(P)1 2 . See Corollary [2] in section 13.21 For 
polygons in § 2 we have Gp(Z) = c + (P)l + d + (P)f(l) where / is a uni- 
versal periodic function. See Corollary [3] in section 13.41 For polygons 
in H 2 we have Gp(l) = c_(P) cosh/. See Corollary H] in section [331 
The coefficients in these formulas depend on how many corners P has 
and on the number of obstacles in its interior. 

Section HJ again, is quite general. In this section we obtain relation- 
ships between partial complexities with one-dimensional base sets and 
counting functions. The main result of this section is Proposition 
It says that if the bases are one-dimensional, then the difference be- 
tween the partial complexity and the counting function is bounded, as 
time goes to infinity. Other assumptions on the base have to do with 
convexity in the phase space. The framework of this section is that of 
piecewise convex transformations [7]. 

In section [5] we specialize again to polygonal billiards. Combining the 
material of preceding sections, we obtain bounds on the position and 
direction complexities for the billard flow and the billiard map. Here is 
a sample of our results. Let P be a euclidean polygon. Let 6 G S 1 (resp. 
z G P) be any direction (resp. position). Let fdg{n) (resp. h z {l)) be 
the direction complexity for the billiard map (resp. position complexity 
for the billiard flow). Then for almost all directions 9 (resp. for almost 
all positions z) we have fd$(n) = 0(n 1+£ ) (resp. h z {l) = 0(/ 2+e )), 
where e > is arbitrary. See Corollary [6] and Corollary [HJ Let now P 
be a spherical polygon, and let h z (l) be the position complexity for the 
billiard flow in P. Then for almost every z G P there is a C = C(z) 
and arbitrarily large I such that h z (l) < CI. See Corollary [TU1 For 
any e > and almost every z G P we have h z {l) = 0{l 1+e ). See 
Corollary [TTJ 

In the study of polygonal billiards the device of unfolding billiard 
orbits is indispensable [I]. If P C M, and (3 is a billiard orbit in P, 
its unfolding is a geodesic in M. Several arguments in section [3] use 
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the technique of lifting billiard orbits to the universal covering space 
of This notion was not written up in the billiard literature. In our 
Appendix section [6] we present the relevant definitions and propositions. 
Proposition [6] puts forward the main property of the universal covering 
space of a geodesic polygon. It relates the unfoldings and the liftings of 
billiard orbits. The proofs in section[3]use Corollary [T3lof Proposition^ 
which deals with the pullbacks of lebesgue measures under unfoldings. 

In order to put our results into perspective, we will now briefly sur- 
vey the literature on billiard complexities. The subexponential growth 
of (full) billiard complexity for arbitrary euclidean polygons is estab- 
lished in [3] and [B]. Both proofs are indirect, in that they do not yield 
explicit subexponential bounds. On the other hand, for rational eu- 
clidean polygons the complexity is cubic. This is contained in [1] for 
convex and in [7J for all rational polygons. The arguments in [T| and 
[7J rely on a theorem in it says that the number of billiard orbits 
between any pair of corners in a rational polygon grows quadratically 
in length. From our viewpoint, this is a statement about the position 
counting functions g z {l). It says that g z {l) = 0{l 2 ) if P C 1R 2 is ra- 
tional and z G P is a corner. By comparison, our Corollary [S] and 
Proposition O yield that g z (l) = 0{l 2+e ) for any e > and almost all 
z G P where P C M 2 is an arbitrary polygon. The directional com- 
plexity fdg(n) has been studied in [8] and [9]. The work [9] concerns 
the directional complexity for the billiard in a rational, planar polygon 
P. Assume that P is convex. Then [9] derives an explicit formula for 
fdg(n), valid for minimal directions 9. (The set of nonminimal direc- 
tions is countable.) By this formula, fdg(n) = 0(n). On the other 
hand, [5] shows that fdg(n) = 0(n d ) for any PcK 2 and an arbitrary 
9. The degree d in the bound does not depend on 9. Our Corollary [BJ 
estimates the complexity fdg(n) for an arbitrary polygon PcK 2 . It 
says that fdg(n) = 0(n 1+£ ) for any e > and almost all directions 9. 

It is plausible that the bounds like Corollary [HI Corollary [HI etc hold 
for any point zGP, any direction 9 G S 1 , etc. 

1. Averages of counting functions 

In this section we introduce the framework of counting functions 
in differentiable dynamics. We will apply it to the billiard dynamics 
later on. Our phase spaces are "manifolds". By this we will mean 
compact manifolds with boundaries, corners, and singular points, in 
general. Our setting involves i) a foliation of the phase space by closed 
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submanifolds that are fibers for a projection onto a manifold of smaller 
dimension; ii) a submanifold in the phase space, transversal to the 
fibers; iii) a weight function on the product of the phase space and 
the time. See the details below. The dynamics in question may be 
discrete or continuous. We will expose the two cases separately. The 
two subsections that follow are parallel, and the treatments differ in 
technical details. 

1.1. Discrete dynamics. Let T : X — > X, T -1 : X — > X be piecewise 
diffeomorphisms with the following data. 

1. There is a fibration rj : X — > O whose base is a compact manifold 
and whose fibers Rg = f]~ l (9) are compact submanifolds, such that 
dim(Rg) = dim(X) — dim(O). We will use the notation X = Ug^Rg. 

2. There is a closed submanifold, Y C X, dim(F) = dim(0), such that 
for k G -P@ the manifolds T k {Y) are transversal to the fibers Rg. 

3. There is a weight function, i. e., a continuous, non-negative function 
w(x,t) on X x N. The function w may depend only on time, e. g., 
w = Xni the indicator function of [0, n — 1]. 

Remark 1. Condition 2 may be weakened, as follows. 
2'. There is a closed submanifold, Y <Z X, and a set ex C O of measure 
zero such that for k G — N and 9 G 0\0 ex the manifolds T k {Y) and Rg 
are transversal. All of our results remain valid if we replace condition 2 
by the weaker condition 2'. However, in our applications to polygonal 
billards, condition 2 may not hold only for polygons in surfaces of 
positive curvature. See section 13.41 To simplify the exposition, we will 
assume in what follows that G ex = 0- 



In view of condition 2, Y{9) = {(x, k) : x G Rg,k G N, T k (x) G Y} is 
a countable (at most) set. The sets Tk(9) = {(x, k) : x G Rg,T k (x) G 
Y} are finite for all keN, and T{9) = UT k (9). 

We define the weighted counting function by 



(x,k)er(e) 

The pure counting function g n {9) corresponds to the weight w = x- 
We have 



(1) 





n 



(2) 




k=0 



2 By convention, N = 0, 1, 
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Proposition 1. Let d6,dy be finite, lebesgue- class measures on 0, Y 
respectively. Then for k G N there are functions > on Y , 

determined by the data 1) and 2) alone, such that 

(3) [ g(6;w)d6 = [ | V>(T- fe ■ y,k)r k (y)\ dy. 
J ® Jy UeN J 

Proof. For any k G N set fk = i] ° T : Y — > O. By conditions 1 and 
2, fk is a local diffeomorphism. Therefore fl{d9) = r k {y) dy. 

It suffices to establish equation ([3]) for the special case w(x,i) = if 
i 7^ k. A point i£l contributes to the integral in the left hand side 
of equation ([3]) iff T k -x G Y, or equivalently, ?7(x) = fk(y), V G F . The 
claim follows by a straightforward change of variables. ■ 

1.2. Continuous dynamics. Let b l : ^ — > ^ be a flow of piecewise 
diffeomorphisms on a phase space \1/ with the following data. 

1. There is a fibration q : — > Z with a compact base and fibers 
q~ 1 (z) = R z C ty, transversal to the flow. We will use the notation 
* = U zeZ R z . 

2. There is a closed submanifold, M C ^ dim(M) = dim(Z) — 1, 
transversal to the flow, and such that iV = U( e iR_6* • M is transversal 
to the fibers i? 2 Jl 

3. There is a weight function, i. e., a continuous, non-negative function 
w(x,t) on \1> x K + . In a special case, u> depends only on time, e. g., 
w = xii the indicator function of [0, 1]. 

In view of condition 2, G(z) = {(x,t) : x G R z ,0 < t,b t (x) G M} is 
a countable (at most) set. The sets Gi(z) = {(x,t) : x G i? 2 ,0 <t< 
I, b\x) G M} are finite for all / G R+, and = \JG t {z). 

We define the weighted counting function by 

(4) g(z;w)= 

(x,t)eG(2) 

The pwre counting function gi{z) corresponds to the weight w = xi- 
We have 

(5) s,(z) = |G,(z)|. 

Proposition 2. Lei <iz, dm &e finite, lebesgue- class measures on Z, M 
respectively; let dt be the lebesgue measure on R. Then there exist a 

3 Our results remain valid if the set of parameters Z cyi C Z where the transversal- 
ity fails has measure zero. See Remark [TJ In what follows, by condition 2' we will 
mean the weakened condition 2 either in the setting of section ITT21 or section [01 
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positive function r(-) on M x R + , determined by the data 1) and 2), 
and such that 

(6) / g(z;w)dz = / • m,t)r(m,t)dmdt 

JZ JMxR + 

Proof. We define the mapping / : M x IR + — ► Z by / = g o ft - *. 
By conditions 1 and 2, / has full rank almost everywhere. The pull- 
back by / of dz is absolutely continuous with respect to dmdt, hence 
f*(dz) — r(m,t)dmdt. 

For < / set Wi(x,t) = w(x,t)xi(t), and let gi(z;w) be the corre- 
sponding counting function. Set hiw) = f z gi(z; w)dz. A point, i£$, 
contributes to h{w) iff x G <p(M x [0, 1]). Under the change of variables 
dz = d(q o <p(m, t)) = r(m, t)dmdt, we have 

Ii(w) = / w(b~ l ■ m,t)r{m,t)dmdt. 

JMx[0,l] 

In the limit / — > oo, we obtain the claim. ■ 

1.3. Special cases. We will discuss a few special cases of Proposition [1] 
and Proposition [3 First, the discrete version. The function g n (9) 
counts the number of visits in Y of points x G Re during the first n 
steps of their journey. Set p k = J Y r k (y)dy, and R n = Ylk=o Pk- Then 
Pk is the volume of Y k = T~ k (Y) with respect to the measure rj*(d9). 
Proposition CD yields 

(7) / g n (9)d9 = R n . 

Je 

In the continuous case the function gi(z) counts the number of visits 
in M of orbits b l ■ x, x G R z , during the period < t < I. Let R(l) be 
the volume of the manifold iVj C ^ with respect to the measure q*(dz). 
Proposition [2] yields 

(8) / gi(z)dz = R(l). 

Jz 

2. Bounds on counting functions 

In this section we analyze the setting of section [T] from the measure 
theoretic viewpoint. This allows us to obtain pointwise upper bounds 
on counting functions in a broad spectrum of situations. 

Let X, ii be a finite measure space. Let f(x;t) (for t G 1R+) be a 
family of nonnegative L 1 functions on X. Set 

(9) Fit) = [ f(x]t)dti(x). 

Jx 



s 



EUGENE GUTKIN AND MICHAL RAMS 



Lemma 1. For almost every x G X there exists C = C(x) > such 
that for arbitrarily large n G N there ist > n satisfying f(x, t) < CF(t). 

Proof. For < C and n G N let 

B n (C) = {xeX: CF(t) < f(x; t) W > n}, 

and set 

B(C) = |J B n {C). 

neN 

Integrating the inequality above, we obtain fi(B n (C)) < C~ l for any n. 
Thus fj, (B(C)) < C~\ and hence /i(nc7 eR+J B(C)) = 0. But n CeR+J B(C) C 
X is the complement of the set of points x G X satisfying the hypoth- 
esis of the lemma. ■ 



Let the setting be as in Lemma [TJ In addition, we suppose that 
i) the functions f(x\ t) are nondecreasing in t and ii) F(t) — > oo. 

Lemma 2. Let e > be arbitrary. Then for almost every iGl there 
exists T = T(x, e) > such that for allt > T we have 

(10) /(x;t)<F(t)(l + log(l+F(t))) 1+£ . 

Proof. Denote by f(x;t~) (resp. F(t~)) the limits of f(x;s) (resp. 
F(s)), as s -> t from the left. For n G N set t n = M{t : F(t) > 2 n }. 
Then F(t( n+1 )~) < 2F(t n ). Let y4 n C X be the set of points satisfying 
the inequality 



(11) f(x;t-) < ^F(t n -)(l + log(l + iF(t n -))) 1+e . 

It suffices to prove that the set UneN flfc>n A-k has full measure. Indeed, 
for x G A n and t G [t n , t n+ i) we have 

f(x;t) < f(x;t- +1 ) < iF(t n+ r)(l + log(l + ^F(t n+ r)) 1+e 

< F(t n )(l + log(l + F(t n ))) 1+£ < F(t)(l + log(l + F(t))) 1+e . 

Thus, the points x G UneN f\>n ^fe have the property equation (1101) . 

If £> n C X is any sequence of sets, we set limsup^^ B n = f] i&N U,>i ^j- 
Let B n be the complement of A n in X. Then lim sup n ^ oc B n is the com- 
plement of UneN f\->n ^-fe- ft remains to prove that //(lim sup n ^ oc £>„) = 
0. 

By Tchebysheff inequality, we have 

(12) f,(B n ) < 2(1 + log(l + iF(t n ")))~( 1+£ ). 
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Set n n = fi(B n ). Suppose first that F is a continuous function. Then 
F(t n ~) = F(t n ) = 2 n . By equation (jT2j) 

/Un <2(l + log(l + 2"))"( 1+£ ), 

hence the series converges. Since 

oo 

//(limsup-Bj < 22 

n—*oo 

no 

for any n G N, the claim follows. 

In general, F need not be continuous. It is thus possible that t n = 
t n+ i for some n G N, implying B n = B n+ i. From the series Yl f-n 
we drop the terms /i n such that B n = B n _\. By equation (fT2l) . the 
remaining terms satisfy 

yU „<2(l + log(l + 2"" 2 ))^ 1+£ ). 

Now the preceding argument applies. ■ 

In sections El [5] we will apply these results in the billiard setting. In 
section [3] we will estimate the integrals equation ([9]), hence the bounds 
provided by Lemmas [U [2] will be more specific. The propositions below 
anticipate these applications. 

Proposition 3. Let the setting and the assumptions be as in Lemma 
[H Let < e be arbitrary. 

1. Let F(t) = 0(t p ) for < p. Then for almost every x G X we have 
f(x;t) = 0(tP+*). 

2. Let F(t) = 0(e at ) for < a. Then for almost every x G X we have 
f{x;t) = 0(e( a+£ )*). 

Proof. The first claim is immediate from Lemma [2] and (logt) 1+e = 
o(t £ ). The second claim follows the same way from t 1+£ = o(e et ). m 

For applications to the billiard map we need a counterpart of Proposi- 
tion[3]for integer-valued time. We state it below. Its proof is analogous 
to the proof of Proposition [3j Moreover, the discrete time case may be 
formally reduced to the continuous time case. We leave details to the 
reader. 

Let X, n be a finite measure space. Let f(x; n), n G N be a sequence 
of nonnegative L 1 functions on X such that for every x G X the nu- 
merical sequence f(x;n) is nondecreasing. Set F{n) = j x f(x;n)dfi. 

Proposition 4. Let < e be arbitrary. Then the following claims 
hold. 



10 



EUGENE GUTKIN AND MICHAL RAMS 



1. Let F{n) = 0{n p ) for < p. Then for almost every x G X we have 
f(x;n) = 0(nP +£ ). 

2. Let F(n) = 0(e an ) for < a. Then for almost every x G X we have 



Remark 2. All of the bounds /(•) = O(-) in preceding propositions 
are equivalent to the formally stronger bounds /(•) = o(-). 

3. Counting functions for polygonal billiard 

We will now apply the preceding material to the billiard dynamics. 
Our billiard table will be a geodesic polygon either in the euclidean 
plane M. 2 , or the round sphere § 2 , or the hyperbolic plane H 2 . We refer 
to [4] , [7] , and section [6] for the background. 

3.1. Direction counting functions for billiard maps in euclidean 
polygons. Let P C I 2 be a euclidean polygon, and let T : X(P) — > 
X(P) be the billiard map. Elements of the phase space X = X(P) 
are oriented geodesic segments in M 2 with endpoints in dP. A segment 
x G X ending in a corner of P is singular, the element Tx is not well 
defined. A billiard orbit x, Tx, . . . , T k ~ l x is a singular orbit of length 
k if T k ~ x x is the first singular element in the sequence. 

Assigning to x G X its direction, r\{x) G S 1 , we obtain a fibration 
f] : X — > S 1 with fibers Rg C X. See figured! We define the counting 
function gd$(n) for singular orbits in direction 9 as the number of phase 
points x G Re that yield singular orbits of length k < n. 

Theorem 1. Let P C M 2 be an arbitrary polygon. Let K(P) be the 
set of its corners. Let a{v) be the angle of v G K(P). Let d6 be the 
lebesgue measure on S 1 . 



Proof. It suffices to prove the claim for a singleton, K = {v}. Let 
Y = Y(v) C X be the set of segments x G X ending at v. Let dy be 
the angular measure on Y. These data fit into the setting of section [TTT] 
and gdein; v) is the pure counting function. 

Let B(z, a) be a conical beam of light with apex angle a emanating 
from z G M 2 . After reflecting in dP, it splits into a finite number of 
beams B(zi,ai) satisfying Yl a i = a - The preservation of light volume 
is due to the flatness of dP. 



f(x;n) = 0{e^ a+£ >). 



Let K C K(P). Then 



(13) 
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By preceding remark, the functions Tk(-) of Proposition [T] satisfy 
rfc(-) = 1. The claim now follows from the special case of Proposition [T] 
considered in section 11.31 ■ 

Let p, q be the numbers of corners, obstacles in P. Let k(P) = 
p + 2q — 2. Thus, P is simply connected iff q = iff k{P ) = p — 2. 

Corollary 1. Let Pel 2 be an arbitrary polygon. Then 

(14) / gd e (n)d6 = Tin(P)n. 

Js 1 

Proof. Follows from Theorem [1] via Y2veK(P) a ( v ) = + 2g — 2)n. m 




Figure 1. Base sets for billiard counting functions 
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3.2. Position counting functions for billiard flows in euclidean 
polygons. Let P C M? be a polygon, and let 6* : — > \& be the 
billiard flow. See section [6] for details. For z G P and v G K(P) let 
<7C Z (Z; be the number of billiard flow orbits that start from z G P and 
wind up at v by time I. Then gc z (7) = J2 veK r P -\ gc z (l; v) is the number 
of singular billiard orbits of length at most / starting from z. This is 
the position counting function for the billiard flow in P. 

Theorem 2. Let P C M 2 be a euclidean polygon, and let dz be the 
lebesgue measure on P. Then for any K C K(P) we have 



Proof. It suffices to prove the claim for K = {v}. We view elements 
of \P as pairs z, 9 where z G P is the basepoint, and 9 is the direction. 
Let M = {(v,9) : (v, -9) G \&}. Let g : * -> P be the obvious 
projection. Its fibers R z are the base sets for the counting functions 
gc z (l;v). See figured! Set w = xi- These data satisfy the assumptions 
of Proposition [21 and gc z (l;v) is the pure counting function. 

We set dm to be the angular measure, and compute the function 
r(m,t) in equation (JSJ). By Corollary [TBI in section [6J r = txi- Propo- 
sition [2] implies the claim. ■ 

When K = K(P), the left hand side in equation ( fT5j) is the average 
of the position counting function. The argument of Corollary [TJ yields 
the following. 

Corollary 2. Let Pel 2 be an arbitrary polygon. Then 



3.3. Position counting functions for billiard maps in euclidean 
polygons. We will now discuss two billiard map analogs of the pre- 
ceding example. Let P C M 2 be a euclidean polygon, and let T : 
X(P) —>■ X(P) be the billiard map. The phase space X = X(P) 
consists of pairs (s, a) where s is the arclentgh parameter on dP, and 
< a < Ti is the outgoing angle. See [U [7] and section [6] for details. 
An orbit x, Tx, . . . , T k ~ l x is singular, of (combinatorial) length k if its 
last segment ends at a corner of P. 

Let s G dP, v G K(P). Define GD s (n;v) to be the set of phase 
points (s, a) G X whose orbits of length less than or equal to n end at 



(15) 




(16) 
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V. Set 



gd s (n;v) 



GD s (n;v)\, god s (n;v) 



E 



sm a. 



(s,a)eGD a (n;v) 



The expressions 
gd s (n) - 



^2 9d s (n;v), god s (n) 

v£K(P) 



^2 god s (n;v) 



v€K(P) 



are the pure position counting function and the optical position counting 
function for the billiard map in P. 

Let z G M 2 and let 7 C R 2 be an oriented piecewise C 1 curve. 
Denote by d z s the projection of the arclength form ds of 7 onto the 
direction perpendicular to the line from z to s G 7. The integral 
J^d z s = I opt (7, z)\ < I7I is the optical length 0/7 viewed from z. 

Let z G P. Unfolding /c-segment billiard orbits emanating from z, 
we obtain a set of linear segments in M 2 . Let d z [P\ k) C M 2 be the 
curve traced by their endpoints. We say that d z (P;k) C M 2 is the 
outer boundary of P, as viewed from z, after k iterates. 

Theorem 3. Let P be a euclidean polygon, and let K C K(P) be a set 
of corners. Then 



Proof. It suffices to prove the claims for a singleton, K = {v}. Let r) : 
X — > dP be the natural projection. Using the arclength parametriza- 
tion, we identify dP with the interval [0, \dP\] Cl For < s < \dP\ 
let R s = rj~ 1 (s) C X be the fiber. Then R s are the base sets for 
the counting functions gd s (n;v), god s (n;v). See figured! Let Y = 
Y(v) C X be the set of phase points whose T~ 1 -orbits emanate from 
v. The assumptions of section 11.11 are satisfied. The weight functions 
are w(s, a, t) = Xn(t) and w (s, a, t) = sin a ■ Xn(t) for the two cases at 
hand. Let <p be the angle parameter on Y. The measures on dP and 
Y have densities ds and dip respectively. 

The integrals in the right hand side of equation ([3]) are over the curves 
d v (P; k), < k < n — 1. The integrands are ds((p) and sin a • ds(tp) = 
d v s(ip) in respective cases. ■ 

We will need estimates on lengths and optical lengths. 



(17) 




(18) 
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Lemma 3. For any polygon PcK 2 there exist < c\ < C2 < oo such 

that for n sufficiently large 

(19) 

n n 

cm 2 < £I>|rf(W*))l < c 2 n\ Cl n 2 < ^2^2m(P,k))\. 

veK k=l vEK k=l 

Proof. There exist positive constants di,d 2 and niQ e N, such that 
for any orbit 7 of the billiard map with m > segments, we have 

d\Y)\ < m < c?2 It I P9- 

Let v G K(P). We will estimate J2k=m \°Pt(d v (P, k))\, as n —> 00. 
Let 81 < 9 < 9 2 be the angular parameter for orbits emanating from v; 
let r(9) be the geometric length of the orbit. Suppose that r\ < r{9) < 
r%. Then the optical length in question is sandwiched between the 
lengths of circular arcs of radii ri, r 2 of angular size 9 2 —9\. By preceding 
remarks, if k is sufficiently large, the bounds r\,r 2 are proportional to 
k. The total angular size does not depend on k. Hence, for sufficiently 
large k we have linear upper and lower bounds on J2 v gk \°pt(d v (P, k))\. 
The other inequality follows from |opt(<9„(P, k))\ < \d v (P, k)\. m 

3.4. Position counting functions for billiard flows in spherical 
polygons. The study is analogous to the planar case discussed in sec- 
tion I3.2( we will use the same notation whenever this does not lead to 
confusion. We denote by dz the lebesgue measure on S 2 , and by a(v) 
the angle of a corner of P. Set 

2 

(20) £(x) = 1 — cos 2 x. 

Theorem 4. Let P C S 2 be a geodesic polygon, and let K C K(P). 
Then 

(21) f Y,9c z (t,v)dz= ($>(")) (-1 + C{l-*\IM) 

Proof. It suffices to prove the claim when K = {v}. Let M = M(v) C 
^ be as in section I3.2[ and let da be the angular measure on it. The 
assumptions 1, 3 of section [1] are satisfied; the transversality of 6* • M 
and R z may fail for at most a countable set of parameters P ex C P. See 
Remark [3] in section [61 Hence, condition 2' is fullfilled, and the results 
of section [L~2l hold. The function gc z (l; v) is a pure counting function. 
The claim now follows from Proposition [2] and Corollary [03 ■ 



Let k(P) be as in section 13.11 
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Corollary 3. Let PcS 2 be an arbitrary polygon. Then 
(22) J gc z (l)dz = (k(P)tt + area(P)) 0/ + ((I - n ) 

Proof. For a spherical polygon we have Yl v eK(p) a ( v ) = area (-P) + 
k(P)tt. Substitute this into equation fl2Tj) . ■ 

3.5. Position counting functions for billiard flows in hyper- 
bolic polygons. Our treatment and our notation are modelled on 
section I3~4l We denote by dz the lebesgue measure on H 2 , and by a(v) 
the angles of corners. 

Theorem 5. Let P C H 2 be a geodesic polygon, and let K C K(P). 
Then 



(23) / gc z (l; v )dz = J a(v) J cosh/. 

J P veK \v€K J 



Proof. We repeat verbatim the proof of Theorem HI and use claim 2 in 
Corollary [ 



Let k(P) be as in section 13.11 
Corollary 4. Let Pel 2 be a polygon. Then 

(24) J gc z {l)dz = (k(P)tt - area(P)) cosh/. 

Proof. Repeat the argument of Corollary [3j use the formula ^ veK r P \ a{v) 
K,(P)7i — area(P) relating the angles and the area of geodesic polygons 
in H 2 . ■ 

4. Relating partial complexities and counting functions 

In this section we establish a framework that will allow us to study 
the complexity of a wide class of dynamical systems. Our motivation 
comes from the billiard dynamics. In fact, polygonal billiard is the 
target of applications for our results. The framework is more general, 
however. The following observations served as our guiding principles. 
First, natural partitions of the billiard-type systems are geared to the 
singularities. Second, the billiard dynamics satisfies a certain convexity 
property that is instrumental in the study of complexity. These princi- 
ples are manifest in the framework of piecewise convex transformations 

There are two approaches to the billiard dynamics: The billiard flow 
and the billiard map. See section [6j The framework of piecewise convex 
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transformations is geared to the billiard map. We begin by establishing 
its counterpart for flows. 

4.1. Piecewise convex transformations and piecewise convex 
flows. A piecewise convex transformation is a triple (X, T, T), where X 
is a two-dimensional convex cell complex, r C X is the graph formed by 
the union of one-cells, and T : X — > X is an invertible map, regular on 
the two-cells of the complex, and compatible with the convex structure 

Let \& be a compact manifold, with boundary and corners, in general. 
Let 6* : \& — * \& be a flow, possibly with singularities; let X C ^ be 
a cross-section. We will assume that the singular set of the flow is 
contained in X. For z G X let t+(z),T-(z) be the times when z G \& 
first reaches X under 6*, 6~* for < t. We assume that for any z G 
there is < e = e(z) such that 6* (2) is regular for |i| < e. 

A piecewise convex flow is determined by the following data: A flow, 
6* : — > ty, a cross-section, X C \&, and the structure of a convex cell 
complex on X, compatible with the poincare map. Billiard flows for 
polygons on surfaces of constant curvature are piecewise convex flows 

m- 

4.2. Partial complexities for maps and flows. Let (A,r n ,T n ) be 
the iterates of a piecewise convex transformation (X, T, T)Jj Let F(T n ) 
be the finite set of open faces of T n ; these are the continuity regions for 
T n . The function f(n) = \F(T n )\ is the (full) complexity of (X,T,T). 

Let R C X be a closed subset. Set 

F R (n) = {AeF(r n ):AnR^<l>}. 

Definition 1. The function = \Fn(n)\ is the partial complexity 

of the piecewise convex transformation (X,T,T) based on the subset 
R. 

Let b l : ^> — * \& be a piecewise convex flow, and let R C ^> be a closed, 
convex set transversal to the flow. For < I let Or{1) be the set of 
regular flow orbits of length / starting from R. Let ao,ai G Or(1). A 
homotopy is a continuous family of regular orbits a p G Or(1), < p < 
1, interpolating between ao, cm- We will say, for brevity, that the orbits 
ao, ai are R-homotopic. We denote by Hr{1) the set of i?-homotopy 
classes. 

Definition 2. T/ie function Hr{1) = \Hf>(l)\ is the partial complexity 
(based on R) of the piecewise convex flow 6* : ^ — > 



They are piecewise convex transformations as well [7]. 
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4.3. Partial complexities and counting functions. In what fol- 
lows we assume that R C is a convex graph without isolated vertices. 
For x G R its valence val(x) is the number of edges of x minus one. 
In particular, if x is an interior point of an edge, then val(x) = 1. Set 
val(i?) = max x€ R val(x). We endow R \ {x} with the graph structure 
where x is replaced by 1 + val(x) vertices; each of them is the endpoint 
of a unique edge. If x, y, z, . . . G R are distinct points, then the induc- 
tively defined graph structure on R without x,y, z, . . . does not depend 
on the order of removing these points. We will denote this graph by 
R \ {x, y, z, . . .}. See figure [2] for an illustration. 

Let E(R) and V(R) be the sets of edges and vertices, and let c(R) 
be the number of connected components of the graph. Let hi = hi(R) 
be the betti numbers of R, and set x(R) — ~ \E{R)\- Then 

c(R) = h Q , x(R) = h - hi. 

Lemma 4. Let R be a finite graph, and let Xi, . . . ,x p G R be distinct 
points. Then 

v p 
(25) x{R) + Yl val ( x i) < c(R \ {xi, . . . , x p }) < c{R) + Y va K x i)- 

i=l i=l 

If R is a forest, then the bound on the right in equation (|25|) becomes 
an equality. 
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Proof. It suffices to prove the claims when R is connected, and we 
remove a single vertex, x. Equation (|25p becomes 

(26) x(R) + val(x) < c(R \ {x}) < c(R) + val(x). 

We have \V(R \ {x})\ = \V(R)\ + val(x), \E(R \ {x})\ = \E(R)\, and 
x(R\{x}) = x(-R)+val(x). Equivalently, we have x(R\{ x }) — h (R) + 
val(x)-hi(R) and h (R\{x}) = h {R)+val{x)+{h 1 {R \ {x}) - h^R)). 
The former (resp. latter) identity implies the left (resp. right) inequal- 
ity in equation ( l26i) . 

When R is a tree, we have c(R \ {x}) = c(R) + val(x), and the 
remaining claim follows. ■ 

We will introduce counting functions for singular orbits of the billiard 
map and the billiard flow. 

By definition, an orbit a = {6*(z), < t < /}, does not pass through 
singular points in It is regular if it does not contain any singu- 
lar points in it is singular if one of its endpoints is singular. The 
set Sr{1) of singular orbits of length at most /, based in R, is finite. 
The quantities gc R {l) = \Sr(1)\ and gdji(n) = \RD T n \ are the count- 
ing functions for singular orbits based in R for the flow and the map 
respectively. 

Now we will relate partial complexities and counting functions. We 
do this for a piecewise convex flow b l : \1/ — > ^> and for a piecewise 
convex transformation (X, T, T) . In both cases the partial complexity is 
based on a 1-dimensional subset, say R. Recall that gcn(l), gdn(n) are 
the respective counting functions, and hn(l), fn(n) are the respective 
complexities. We will refer to these situations as the continuous case 
and the discrete case respectively. 

Proposition 5. Let the setting be as above. Then the following state- 
ments hold. 

1. In the continuous case there exist h £ N and l £ IR + such that 
hail) — ho + genii) for l < I. 2. In the discrete case there exist 
fo, no £ N such that for n < n we have fn(n) = fo + gd^n). 

Proof. In both cases the graph R is equipped with a tower of finite sets, 
say X(l) and X n respectively. Let Xoo C R be their union. We will 
compare the number of connected components of graphs R\X(l), R\X n 
with the cardinalities of these sets. 

We consider the discrete case, leaving the continuous case to the 
reader. Let m < n be any pair of natural numbers. By (the proof of) 
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Lemma HI 

c(R\X n )-c(R\X m ) = [h 1 (R\X n )-h 1 (R\X m )]+ £ val(z). 

f"f?. 

We have \ X n ) < h±(R\ X m ); the inequality holds iff X n \ 

X m breaks cycles in R \ X m . Since the sequence hi(R \ X k ) G N is 
nonincreasing, it stabilizes. Thus, there exists n\ G N such that for 
rii < m < n we have hi(R \ X n ) = h 1 (R\ X m ). 

The set of points i6i2 satisfying 1 < val(x) is finite. Thus, there 
exists ?7,2 G N such that if ri2 < k and x G X^ \ Xk, then val(x) = 1. 

Set no = max(ni,ri2). Then for < m < n the above equation 
yields c(R \ X n ) — c(R \ X m ) = \X n \ X m \. Specializing to m = n , we 
obtain f R (n) = (f R (n ) - gd R (n )) + gd R (n). m 

5. Bounds on partial complexities for the billiard 

We will use the preceding material to derive bounds on partial com- 
plexities for the polygonal billiard. 

5.1. Direction complexities for billiard maps in euclidean poly- 
gons. We use the setting and the notation of section 13. 11 For a polygon 
P and a direction 8, we denote by fdo(n) the partial complexity with 
base Re. This is the complexity in direction 9. 

Corollary 5. For lebesgue almost all directions 9 there is C = C(9) 
and there are arbitrarily large n such that fdg(n) < Cn. 

Proof. Each Rg is a convex graph in the phase space [7] . By Lemma Q] 
and Corollary [U the counting functions gdg(n) have the desired proper- 
ties. By the second claim of Proposition the directional complexities 
do as well. ■ 

Corollary 6. For any e > and almost every direction 9 we have 
fd e {n)=0{n 1 ^). 

Proof. The proof goes along the lines of the proof of Corollary In- 
stead of Lemma [U we use Proposition |4] (the first claim). ■ 

5.2. Position complexities for billiard flows in euclidean poly- 
gons. Let P be a euclidean polygon, and let z G P be any point. We 
consider the billiard flow in P, and use the setting of section [3721 Thus, 
gc z {l) is the position counting function for orbits emanating from z. We 
denote by h z {l) the corresponding partial complexity. 

Corollary 7. For almost every point z there is a positive number C = 
C(z) such that h z (l) < CI 2 for arbitrarily large I. 
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Proof. The sets R z satisfy the assumptions of section HI The claim 
follows from Lemma dj Corollary [2] and the continuous case in Propo- 
sition [51 ■ 

Corollary 8. For any e > and almost every z G P we have h z {l) = 
0(l 2+e ). ' 



Proof. The proof is similar to the preceding argument, and we use the 
first claim in Proposition [3] instead of Lemma [TJ ■ 

5.3. Position complexities for billiard maps in euclidean poly- 
gons. This is the billiard map analog of the preceding example. Let 
P be a euclidean polygon, and let s G dP. We use the setting of sec- 
tion [331 There we have defined the counting functions gd s (n), god s (n). 
Let f s (n) be the partial complexity corresponding to gd s (n). This is 
the position complexity for the billiard map. 

Corollary 9. Let P C M 2 be a polygon such that ^2 V&K Ylk=i \^v(P] k)\ 
has a quadratic upper bound^ Then for almost all s G dP we have 
f s (n) = 0{n 2+£ ) for any < e. 

Proof. The sets R s C X satisfy the assumptions of section [U We use 
Theorem [31 Lemma [31 and apply Proposition [51 ■ 

The estimate of Corollary [9] on f s (n) is conditional, because in gen- 
eral we have no efficient upper bound on Ylk=i \dv(P] k)\. 



5.4. Position complexities for billiard flows in spherical poly- 
gons. We use the setting of section 13.41 For a spherical polygon, 
P C S 2 , and z G P, let h z (l) be the position complexity. 

Corollary 10. For almost every point z G P there is C = C(z) and 
there are arbitrarily large I such that h z (l) < CI. 

Proof. The sets R z satisfy the assumptions of section[H We use Lemma[U 
Corollary [31 and Proposition [51 ■ 

Corollary 11. For any e > and almost every z G P we have h z [l) = 
0(/ 1+£ ). ' 

Proof. See the proof of Corollary [HI ■ 



5 This holds if P is a rational polygon [TT] . 
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5.5. Position complexities for billiard flows in hyperbolic poly- 
gons. This material is the hyperbolic plane counterpart of section 
and we use the setting of section 13.51 

Corollary 12. Let P C H 2 be a geodesic polygon, let z G P , and let 

h z (l) be the position complexity. Then for almost every point z G P we 
have h z {l) = 0(e( 1+£ ) / )- 

Proof. We verify that the sets R z satisfy the assumptions of section [4j 
and mimick the proof of Corollary [HJ we use Corollary HJ Proposition [31 
and the continuous case of Proposition [5j ■ 

6. Appendix: Covering spaces for polygonal billiards 

Let M be a simply connected surface of constant curvature, and let 
P C M be a connected geodesic polygon. We normalize the metric so 
that the curvature is either zero (M = M 2 ), or one (M = § 2 ), or minus 
one (M = M 2 ). 

Let A be the set of sides in P. We will denote its elements by a, b, ... . 
For a side, say a G A, let s a G Iso(M) be the corresponding geodesic 
reflection. We associate with P a Coxeter system (G, A) [2]. We denote 
by a a , (Tfe, • • • G G the elements corresponding to a, b, ... G A. They 
generate G. The defining relations are a\ = 1 and (o- a o~b) n( ' a ' b ^ — 1; 
the latter arise only for the sides a, b with a common corner if the 
angle, 9(a,b), between them is 7r-rational. In this case n(a,b) is the 
denominator of 9(a, b)/ir. Otherwise n(a, b) = oo. 

To any "generalized polyhedron" P corresponds a topological space 
£ endowed with several structures, and a Coxeter system [2]. Our 
situation fits into the framework of [2], and we apply its results. First, 
£ is a differentiable surface. Second, £ is tiled by subsets P g ,g G G, 
labelled by elements of the Coxeter group G; we call them the tiles, and 
identify P e with P. The group G acts on £ properly discontinuously, 
preserving the tiling: g ■ P h = P gh . 

Since P e is identified with P C M, it inherits from M a riemann- 
ian structure. The action of G is compatible with this structure, and 
extends it to all of £. This riemannian structure generally has cone 
singularities at vertices of the tiling <t = Ug^cPg^ Around other points 
this riemannian structure is isometric to that of M; in particular, ex- 
cept for cone points, £ has constant curvature. The group G acts on <£ 
by isometries. 



6 Each vertex, v, corresponds to a corner of P. The metric at v is regular iff the 
corner angle is w/n, n = 2, 3, . . . . 
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Definition 3. The space £ endowed with the riemannian structure, 
the isometric action of G and the G-invariant tiling € = Ug^Pg is the 
universal covering space of the geodesic polygon P C M. 

If X is a riemannian manifold (with boundary and singularities, 
in general), we denote by TX = U X&X T X X its unit tangent bundle. 
The classical construct of geodesic flow, G x : TX — > TX, extends to 
manifolds with boundaries and singularities. In particular, G x makes 
sense when X = M, P, or <£. Another classical construct, the expo- 
nential map, also extends to our situation. For i6las above, and 
(v,t) G T X X x R + , we set exp x (v,t) G X be the base-point of G x (v). 
We will use the notation exp^ to indicate that we are exponentiating 
from the point x. If X is nonsingular, then exp^ : T X X x M + — > X is 
a differentiable mapping. For X with singularities, such as our P and 
C, the maps exp^- are defined on proper subsets of T X X x M + ; these 
subsets have full lebesgue measure. Generally, the maps do not extend 
by continuity to all of T X X x M.+. 

Let X, Y be nonsingular riemannian manifolds of the same dimen- 
sion; let if : X — > Y be a local isometry. It induces a local dif- 
feomorphism $ : TX — > TY commuting with the geodesic flows: 
$ o G x = G Y ° The exponential maps commute as well: (f o exp^ = 
expy ^ o d x ip. These relationships hold, in particular, for coverings of 
nonsingular riemannian manifolds. Suitably interpreted, they extend 
to (branched) coverings of riemannian manifolds with boundaries, cor- 
ners, and singularities. In our case X = £, while Y = M, or Y = P. 
We will now define the mappings f : <£ —> P, F : T<£ — > TP and 
r - : IT — .U. «I> : VlT — TM. 

The identification P e = P defines /, tp on P e . To extend them to all 
of £, we use the tiling <£ = U 96 g-P 5 and the actions of G on £ and M. 
In order to distinguish between these actions, we will denote them by 
g ■ x and g(x) respectively. Then there is a unique x G P e such that 
z = g ■ x. We set f(z) = x G P and <p(z) = g(x) G M. By basic 
properties of Coxeter groups [2] , the mappings /, <p are well defined. 
Moreover, / : C — > P and tp : <£ — > M are the unique G-equivariant 
mappings which are identical on P e By construction, both mappings 
are continuous; they are diffeomorphisms in the interior of each tile, 
P g C C, and on the interior of the union of any pair of adjacent tiles. 

The potential locus of non-differentiability for both / and tp is the 
set V of vertices in the tiling £ = U g< z G P g . We have V = f~ l (K(P)) 
where K(P) is the set of corners of P. By equivariance, <p(V) = 



7 The action of G on P is trivial. 
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\J ge ag{K{P)) C M@ There are two kinds of points in V: vertices 
coming from the corners of P with 7r-rational and 7r-irrational angles. 
Their cone angles are integer multiples of 2ir and are infinite respec- 
tively. Vertices v G V with cone angle 2ir are, in fact, regular points 
in £, and the mappings /, <p are both regular there. Around a vertex 
v with cone angle 2k-K > 2n the mapping p is different iable, but not a 
diffeomorphism; it is locally conjugate to z t— > z k . Near such a vertex, 
tp is a branched covering of degree k. At a vertex with infinite cone 
angle, tp has infinite branching. 

Remark 3. The set <p(V) C M is countable. (It is finite iff the 
group generated by geodesic reflections in the sides of P is a finite 
Coxeter group. Typically, tp(V) C M is a dense, countable set.) 
Let M = S 2 , and let 2 i— > 2' denote the antipodal map. Set F = 
P fl (ip{V) U (<y2(V))'). Points of F are exceptional, in the following 
sense. Let z G P be such that the beam P^ of billiard orbits emanat- 
ing from z contains a sub-beam focusing at a corner of P. Then z E F. 
This follows from Proposition below. 

Thus, P contains all points z G P for which the transversality as- 
sumption in Condition 2 of section 11.21 fails. Since F is countable, the 
set of exceptional parameters has measure zero, and Condition 2' is 
satisfied. See Remark [1] in section [TJ 

Furthermore, the mappings / and ip are local isometries. They are 
isometries on every tile P g C <£; we have f(P g ) = P, <p(P g ) = g(P) C 
M. Let p • a be a side of P g , let ft. = a a g and let be the adjacent 
tile. The maps / : P g -> P, P h -> P and ^ : P g -> ^(P), P h -> ft(P) are 
coherent around the common (open) side (? • a. The map / is never an 
isometry on P g UPh', for tp this is the case iff the interiors of g(P), h(P) 
are disjoint in M. The latter generally fails for nonconvex P. 

By coherence of / and p across the sides separating adjacent tiles, 
we lift them to the tangent bundles, obtaining the mappings of unit 
tangent bundles F : T(£ — > TP, $ : T€ — > TM, which are also defined 
on vectors based at the vertices of the tiling £ = U ge cPg- Let v be a 
vertex, and let a be the angle of the corner f(v) G K(P). Then $ : 
T v £ —>■ T,pf v \M is m-to-1 if a = vrm jn and oo-to-1 if a is 7r-irrational. 
The geodesies 7(t) in (t cannot be further extended (generally) once 
they reach a vertex. All other geodesies in <t are defined for — 00 < t < 
00. 



The representation M = U g ^ag{P) is not a tiling, in general. 
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Using the inclusion P C M, we identify TP with the subset of TM 
consisting of M-tangent vectors with base-points in P, and directed in- 
ward. Any v G TP defines the billiard orbit in P, (3{t) = exp P (tv), < 
t, and the geodesic in M, j(t) = exp M (tv),0 < t. They are related 
by the canonical unfolding of billiard orbits. This is an inductive pro- 
cedure which replaces the consecutive reflections about the sides of P 
by consecutive reflections of the "latest billiard table" g(P) about the 
appropriate side, yielding the next billiard table h(P), and continuing 
the geodesic straight across the common side of g(P) and h(P). See jl] 
in the planar case and [7], section 3.1, in the general case. Let i6P 
and let v G T X P. We denote by (3 V (resp. j v ) the billiard orbit in P 
(resp. the geodesic in M) that emanates from x in the direction v. The 
unfolding operator, U x : fi v i— ► j v , preserves the parametrisations. 

Proposition 6. Let x G P, v G T X P. Identify P and P e C £ and let 

x G P e , v G T x € be the corresponding data. Then for t G K + we have 

(27) U x (exp P (v, t)) = (/^expg-O, t)). 

Proof. We will freely use the preceding discussion. As t G M + goes to 
infinity, exp P (v, t) runs with the unit speed along a billiard orbit in P. 
The curve exp £ (w, t) is the geodesic in £ defined by the data (x, v), and 
(p(exp<r(v,t)) is the geodesic in M emanating from x in the direction 
v. The billiard orbit in P and the geodesic in M are related by the 
unfolding operator. ■ 

For x G P let E X P = T X P x IR + be the full tangent space (or the full 
tangent cone) at x. If S C T X P is a segment, let ES X P = S x R + be 
the corresponding subcone. We use the analogous notation for i 6 £ 
or x G M. In polar coordinates (£, 9) in M 2 the lebesgue measure on 
E X P is given by the density tdtdO. 

Corollary 13. Let x G P C M be arbitrary, and let expf, : E X P — > P 
be the exponential mapping. The pull-back by expp of the lebesgue 
measure on P to E X P is the smooth measure with the density dv{t,9). 

1. When M = R 2 , we have du = tdtdd. 

2. When M = B 2 , we have dv = smb. tdtdO. 

3. When M = S 2 , we have dv = \ sint\dtd9. 

Proof. By Proposition the measure in question coincides with the 
pullback to the tangent space E X M of the riemannian measure on M 
by the exponential map E X M — > M. The latter is well known. ■ 

We point out that the preceding material has a billiard map version. 
We will briefly discuss it now. Let = (z(t),9(t)), t G M, be an orbit 
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of the billiard flow. We obtain the corresponding billiard map orbit 
f3d(k), k G Z, by restricting j3(t) to the consecutive times t& such that 
z{tk) G dP. The correspondence /?(•) i— > fld(') is invertible. This allows 
us to formulate the billiard map versions of the universal covering space, 
the lifting of billiard map orbits to the universal covering space, and the 
relationship between the liftings and the unfoldings, a la Proposition [6j 
Since we are not directly using this material in the body of the paper, 
we spare the details. 
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